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ABSTRACT 

Eourier analysis on a cut sky of the two-year COBE DMR 53 and 90 GHz sky maps 
(in galactic coordinates) is used to determine the normalisation of an open inflation model 
based on the cold dark matter scenario. The normalised model is compared to measures 
of large-scale structure in the universe. Although the DMR data alone does not provide 
sufficient discriminative power to prefer a particular value of the mass density parameter, 
the open model appears to be reasonably consistent with observations when Oq ~ 0.3 — 0.4 
and merits further study. 
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1. INTRODUCTION 


Observational evidence (summarized in Ratra & Peebles 1994b, hereafter RPb) mo¬ 
tivates the consideration of a cosmogony with a clustered mass density parameter, Oq, 
significantly less than unity, but possibly somewhat larger than the baryon mass density 
parameter, required for standard nucleosynthesis (e.g.. Reeves 1994). In this Letter^ we 
focus on a low-density cold dark matter (CDM) model with open spatial sections and zero 
cosmological constant A, with an early inflation epoch (Ratra & Peebles 1994a, hereafter 
RPa; RPb). In this model energy-density perturbations are generated by ground-state 
zero-point fluctuations during inflation (Lyth & Stewart 1990; Ratra 1994; RPa).f 

Kamionkowski et ah (1994; hereafter KRSS) previously considered the open model 
discussed here and computed the cosmic microwave background (hereafter CMB) radiation 
temperature anisotropy from the generalized Sachs-Wolfe relation between the radiation 
perturbations and the energy-density inhomogeneity (Gouda, Sugiyama, & Sasaki 1991; 
Sugiyama & Silk 1994). However, KRSS normalised the model to the two-year DMR value 
for the rms temperature anisotropy at 10® angular resolution (Bennett et ah 1994, hereafter 
B94). This KRSS normalisation technique is unsatisfactory since it fails to account for the 
correct DMR beam Alter function and the galaxy cut as applied to the data (Wright et 
ah 1994a) and consequently underestimates the correct normalisation amplitude. The 
low observed quadrupole on the sky also affects the normalisation inferred from the rms 
temperature anisotropy (Banday et ah 1994). 

The appropriate normalisation is determined by using all the angular information 
present in the data, as best represented by the power spectrum, the amplitude of which 
can be characterized by the rms quadrupole anisotropy amplitude, Qrms-PS- Here we 
use the incomplete-sky power spectrum estimation technique (described in Gorski 1994, 
hereafter G94; and applied to the DMR data in the context of a power law GMB anisotropy 
model in Gorski et ah 1994a, hereafter G94a) to determine Qrms-PS, as a function of Oq, 
from the two-year DMR galactic sky maps. With the normalisation fixed, we then proceed 

f In a related model where an open-inflation bubble is created by tunnelling in an 
inflating, spatially-flat, de Sitter spacetime, it might prove more appropriate to use a 
combination of the ground and excited states as an initial condition in the second inflation 
epoch (Bucher, Goldhaber, & Turok 1994); the final spectrum is probably determined by 
the model assumed for the first epoch. 
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to estimate various statistics of cosmological interest and subsequently comment on the 
viability of the open model as a description of our universe. 

2. SUMMARY OF COMPUTATION 

The determination of the CMB power spectrum from the DMR data is hindered by 
the necessity to eliminate data near and on the galactic plane (because of the strong, non- 
cosmological, emission) from the sky maps. The proper Fourier analysis of temperature 
anisotropy on the cut-sky data-set should utilize orthogonal basis functions. Following the 
anisotropy power spectrum technique in G94, a set of orthonormal functions is constructed 
for the 4016 pixels surviving the |6| > 20® galactic cut, and the sky temperature anisotropy 
is decomposed into 961 Fourier coefficients (probing the spectral range I < 30). As in G94a, 
both systematic errors (which are strongly limited, B94) and faint residual high-latitude 
galactic signals are neglected, and the combined 53 and 90 GHz two-year galactic sky 
maps are used in the analysis. The DMR data together with the open model multi¬ 
pole coefficients, parameterized by Qrms-PS and Oq, define an exact gaussian likelihood 
function T((5rms-PS, f^o)- We have computed the GMB anisotropy multipole moments. 
Cl = (|a™P), where 8T/T = ^ for Oq = 0.03,0.06, and in steps of 0.1 from 

Oo = 0.1 to 1.0. The spectra used in this work are shown in Fig. 1, where we have adopted 

= 0.03 and h = 0.5 (the Hubble parameter Hq = lOOh km s“^ Mpc“^). The shape 
of the low-t' CiS is almost entirely governed by Oq and is relatively insensitive to h, 
and ionization history (using the tabulated values of h and = 0.0125/i“^, for t' < 20 
we find \CilC 2 \ is shifted by < 3% for 0.1 < Oq < 0.5 and Oq = 1, and leaves the inferred 
normalisation essentially unchanged). The monopole and dipole are algebraically excluded 
from the analysis (as described in G94). We have analyzed the data both including and 
excluding the observed quadrupole moments, since these are the harmonics most likely to 
be contaminated by residual galactic and systematic emission signals. The results of this 
likelihood analysis are summarized in Figs. 2-4 and Table 1 (columns 3 and 4). Fig. 2 
shows the likelihood surface for the more restrictive case with the 8 = 2 mode included in 
the analysis (with the usual assumption of uninformative, uniform prior in Qrms-PS and 
Oo)- When the quadrupole is excluded the likelihood amplitude shows less suppression 
near Oq ~ 0.15. In both cases the 95% contours of the likelihood extend throughout the 
entire range of Oq and hence no models can be formally rejected. 
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Using the normalisation derived from the DMR data, we determine the rms values 
for a number of measures of cosmological structure (see Table 1). For each value of Oq 
we record in column (2) a value for h which results in an expansion time to 12Gyr 
(the last line is the Einstein-de Sitter case with to ~ ISGyr), consistent with the lower 
end of recent determinations (Renzini 1993; Tonry 1993), and take = 0.0125/i“^, 
the standard nucleosynthesis estimate (e.g.. Reeves 1994). The quoted ranges for the 
tabulated statistics only account for the uncertainty in the inferred model normalisation, 
but ignore the uncertainty in the values of h and as derived above. Such uncertainties 
do not significantly affect the Qrms-PS inference, but do affect estimates of the other 
statistics. Future work should more rigorously account for the allowed parameter range, 
as the accuracy of the Qrms-PS determination improves. For the time being, we normalise 
the model to Qrms-Psf from the no-quadrupole case (Table 1, column 4). 

Golumn (5) gives the present rms linear mass fluctuation, (SM/M), within an 8h~^ 
Mpc sphere, determined from a numerical integration of the perturbation equations, where 
the range here [and in columns (6) - (11)] only accounts for the Icr uncertainty in the esti¬ 
mate of Qrms-PS- Golumns (6), (8), and (10) give the present rms value ^ of \/3 (from 
isotropy) times the (gauge-invariant) radial peculiar velocity (RPa, eq. [11.59]) averaged 
over spheres of radius 40/i“^, Q0h~^ and 100/i“^ Mpc, computed from a numerical integra¬ 
tion of the perturbation equations. However, the observed large-scale velocities (which are 
radial) are conventionally quantified in terms of the rms value Vp of the three-dimensional 
peculiar velocity. Since we cannot express Vp as a sum over modes in an open space, we use 
the flat space expression (Peebles 1993, hereafter PPG, eq. [21.68]), with the appropriate 
open space power spectrum and window function, and give Vp averaged over spheres of 
radius 40/i“^, Q0h~^, and 100/i“^ Mpc in columns (7), (9), and (11). Interestingly, the 
computed values of Vp and Vp are somewhat discrepant, contrary to intuitive expectation 
on the scales considered here (much smaller than the curvature radius). We believe that 
this phenomenon merits further investigation. 

f These numerical values do not significantly alter the reheating redshift estimated in 
RPb with a different assumed Qrms-PS normalisation; we note that unlike some inflation 
models, in this model the observations do not impose a stringent upper bound on the slope 
of the inflation epoch scalar held potential (Ratra 1992; RPa). 
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3. DISCUSSION 


The likelihoods shown in Fig. 3 indicate that the DMR data alone cannot rule out 
any value of Oq (including the quadrupole, Oq ~ 0.15 is acceptable at > 5%, Fig. 2, 
and this is the most restrictive case). While the DMR data is consistent with the low-t' 
CiS in the scale-invariant Einstein-de Sitter model (G94a, Wright et ah 1994b), it is also 
consistent with a flat A-dominated case (Bunn & Sugiyama 1994; Gorski, Stompor, & 
Banday 1995), and, as shown here, with an open model. Thus, the two-year DMR data do 
not necessarily require a flat universe or a scale-invariant model. The diminutive (galaxy- 
subtracted) CMB quadrupole (B94) might be considered evidence against the flat, scale- 
invariant model, and weakly supportive of the open models with suppressed quadrupole 
(Figs. 1 - 3). Nevertheless, it is as yet premature to draw such conclusions, based solely 
on the DMR data, since it is precisely in the very low order multipoles that the cosmic 
variance in all models is large and the galactic foreground contribution remains uncertain. 

The rms fluctuation in the number of galaxies in a randomly placed sphere of radius 
8/i“^Mpc, (SN/N) = 0.79 to 1.1 (PPG, eqs. [7.33, 7.73]). From column (5), with (SN/N) = 
b(SM/M)^ where b is the linear bias factor, Oq $ 0.2 requires an unreasonably high bias, 
Oo ~ 0.3 — 0.5 could be consistent with mild bias, while Oq ^ 0.7 needs antibiasing. 
For Oo ~ 0.3 — 0.5, {SM/M)[8h~^M.pc] more consistently reproduces the observed cluster 
mass and correlation functions (e.g., Bahcall & Gen 1992; White, Efstathiou, & Erenk 
1993; although Oq = 0.3 might require a somewhat higher h and lower Qb than used 
here). Eurthermore, the shape of the observed galaxy power spectrum is better fitted 
by a low-density GDM model. The GDM power spectrum shape parameter corrected 
for fls dependence, (e-g -5 Efstathiou, Sutherland, & Maddox 1990; Peacock & 

Dodds 1994), is determined from miscellaneous galaxy and galaxy cluster catalogs to be 
~ 0.255±0.017. Interestingly, this is consistent with the models in Table 1 for Oq ~ 0.3 —0.4 
and h ~ 0.65 (which is in the range now under discussion — see Ereedman et ah 1994, 
and references therein). 

Observed large scale flows, as currently described in the literature, are problematic 
for the open model. The POTENT analysis of Bertschinger et ah (1990) has rendered 
254 < Up(40/i“^Mpc)/(kms“^) < 522, 163 < Up(60/i“^Mpc)/(kms“^) < 491, (95% c.h). 
Accounting for the cosmic variance in the measurement of one bulk velocity (Bond 1994, 
eqs. [47, 48]) suggests that one need not yet reject any models at Oq ^ 0.3. On scales 
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~ 100/i“^ Mpc the corresponding range is 333 < Up/(kms“^) < 1045 (Lauer & Postman 
1994), implying Oq ^ 0.7. Hence, the observed velocity fields can significantly constrain 
low-density models, however, a more rigorous analysis (that accounts for the survey window 
function and the hypersurface dependence of the velocity) as well as more and better data 
is required to draw a firm conclusion. 

Dynamical estimates of mass clustered on scales ^ 10/i“^Mpc indicate 0.1 $ Oq $ 0.3 
(PPC, §20). The situation on scales ^ 20/i“^Mpc is not yet entirely clear: if galaxies trace 
mass, the IRAS/POTENT analysis indicates Oq > 0.5 at the 95% confidence level (Dekel 
1994), but some other analyses are also consistent with lower Oq (e-g-, Roth 1993; Dekel 
1994). Nusser & Davis (1994) have recently presented an estimation of (= Oo'^ ®/&iras = 
/M)[8h~^M.pc], where 6iras is the linear bias factor for IRAS galaxies) and 
find, at ~95% confidence, 0.2 < /3i < 1.0. This, interestingly enough, corresponds to 
models with 0.3 $ Oq $ 0.6 in Table 1. 

The DMR detection of CMB anisotropy at a level significantly higher than expected 
in the biased scale-invariant fiat CDM model focussed attention on alternative Oq = 1 
models, such as ‘tilted’ CDM or mixed dark matter (MDM) models, which appear consis¬ 
tent with present epoch observations. However, when normalised to the DMR data, these 
models tend to form galaxies too recently. In the tilted case Cen & Ostriker (1993, here¬ 
after CO) find that when {SM/M)[8h~^M.pc] = 0.5 galaxy formation peaks at a redshift 
Zg ~ 0.5 (since CO used a somewhat large (SM/M)^ the upward revision of Qrms-PS will 
not raise Zg significantly), contrary to observational evidence suggesting that the giant 
elliptical luminosity functions at 2 ; = 1 and 0 are not dissimilar (e.g., Cowie, Songaila, & 
Hu 1993). The MDM model appears to predict too few damped Lyman-a systems (e.g., 
Ma & Bertschinger 1994), which are thought to be young galaxies, and are seen out to 
redshifts as high as 2 ; ~ 3 (Lanzetta 1993). It seems encouraging that adiabatic low-density 
models, such as the A-dominated fiat model and the open model considered here, as well 
as low-density isocurvature models (Peebles 1994) and higher-bias low-density topological 
defect models (Pen & Spergel 1994) tend to form structure early. For instance, scaling 
from the CO result for the tilted case with Zg = 0.5, the Oq = 0.4 open model with 
{SM/M)[8h~^M.pc] = 0.8 would have Zg = 2.5 (e.g., RPb), which may be adequate if 
objects like the ‘young cluster’ at 2 ; = 1 (Le Fevre et ah 1994) turn out to be common. 
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In summary, we have used the DMR two-year galactic maps to normalise an open 
inflation model, and have computed large-scale structure statistics in this model. The 
DMR data alone do not strongly constrain Oq, hut somewhat disfavour Oq ~ 0.15 (which 
is consistent with some non-CMB observations). For Oq ~ 0.3 — 0.4 the open model 
appears to be mostly consistent with observations, albeit not too impressively so. More 
detailed modelling (particularly of large scale flows and small angle CMB anisotropies) and 
improved observational constraints will be required to resolve definitively the compatibility 
of this model with the observed properties of the universe. 

We acknowledge the efforts of those contributing to the COBE DMR. We thank D. 
Bond and M. D avis for helpful discussions. This work was supported in part by the Office 
of Space Sciences of NASA Headquarters, by NSF grant PHY89-21378, by the David and 
Lucile Packard Foundation, and by a JSPS Postdoctoral Fellowship for Research Abroad. 
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FIGURE CAPTIONS 

Fig. 1.- CMB anisotropy multipole coefficients for the open inflation model with h = 0.5, 
= 0.03, and Oq = 0.03, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 in ascending 
order, determined numerically from a Boltzmann transfer code. The spectra are normalised 
relative to the C 2 quadrupole amplitude and offset from each other to aid visualisation. 
At low 7 the shape of the Ci spectrum is determined by Oq (which fixes the perturbation 
growth rate and the ratio of the Hubble and space curvature scales), and by the wavelength 
dependence of the energy-density power spectrum. 

Fig. 2.- The likelihood function T((5rms-PS , (normalised, for clarity, to 0.9 at the 
highest peak near Oq = 0.03, Qrms-PS ~ 15/iK) derived from a simultaneous analysis 
of the DMR 53 and 90 GHz two-year data for Fourier modes 7 G [2,30]. Excluding 
the quadrupole results in a qualitatively similar result at large and very small Oq, with 
the likelihood less suppressed at Oq ~ 0.15. The contour plot of the likelihood shows 
the 99.9% (outermost), 95%, and 68% (innermost, disconnected) contour lines, and the 
projected location of the highest likelihood peak. One should note, on comparing to Fig. 
1, that models with suppressed low-7 multipoles have somewhat higher likelihood (at low 
Oo and Oq ~ 0.7) than those with enhanced low-7 power (Oq ~ 0.2). However, the relative 
differences in the amplitude of the likelihood are too small to allow definite rejection of 
any value of Oq- 

Fig. 3.- Likelihood densities for Oq, derived from T((5rms-PS , f^o) (and with arbitrary 
equal-integrated-area normalisation). Light type shows the projected likelihood and heavy 
type the marginal likelihood [oc J <i(5rms-PsR(Qrms-PS , f^o)]. The solid lines are the case 
including the quadrupole, the dotted lines are the quadrupole excluded case. 

Fig. 4.- Gonditional likelihood densities for Qrms-PS, derived from L((5rms-PS, f^o) (and 
with arbitrary equal-integrated-area normalisation), for Oq = 1, 0.5, 0.4, 0.3, 0.2 running 
from left to right. Solid lines are the case with the quadrupole, dotted lines are the no- 
quadrupole case. 
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